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We study the dynamical stability of holographic superconductors. We first classify perturbations 

around black hole background solutions into vector and scalar sectors by means of a 2-dimensional 

rotational symmetry. We prove the stability of the vector sector by explicitly constructing the 

positive definite Hamiltonian. To reveal a mechanism for the stabilization of a superconducting 

phase, we construct a quadratic action for the scalar sector. From the action, we see the stability 

of black holes near a critical point is determined by the equation of motion for a charged scalar 

field. We show the effective mass of the charged scalar field in hairy black holes is always above the 

Breitenlohner-Freedman bound near the critical point due to the backreaction of a gauge field. It 

implies the stability of the superconducting phase. We also argue that the stability continues away 

04 ■ from the critical point. 
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I. INTRODUCTION 



It is believed that the AdS/CFT correspondence is useful to study strongly coupled systems [1(. Remarkably, the 
AdS/CFT correspondence has been extended to the correspondence between gravity and condensed matter physics [2|. 
In particular, it has been shown that there exists a gravity dual of a superconductor |3j-|5|. Interestingly, the mean 
field square root behavior of the Landau-Ginzburg second order phase transition has been reproduced through the 
C"| \ gravity dual of the superconductor. Moreover, it turned out that conductivity has similar features to that of the 
superconductor (see review articles and references therein [y, 0] ) . 

It is well known that the superconductor can be explained as the second order phase transition phenomena. In 

terms of the Landau-Ginzburg theory, below a critical temperature, the effective mass of the order parameter field 

^i ■ becomes tachyonic and destabilize normal phase. Eventually, a superconducting phase is realized as the new phase 

where the stability is guaranteed by a quartic potential. In the holographic description of the superconductor, the 

^^ ■ onset of the instability is well understood. The point is that the coupling of a charged scalar field to a gauge field 

^O | through covariant derivatives induces an effective mass term for the scalar field [8|. This term becomes relevant as 

r~^. ■ temperature gets lowered at fixed charge density, and eventually it makes the effective mass of the charged scalar field 

below the Breitenlohner-Freedman (BF) bound Q, and hence destabilize the system. On the other hand, the fate of 

the instability is not fully understood. In fact, in the gravity dual of superconductors, no quartic potential to stabilize 

the system after the instability is present. Hence, although a trigger of the phase transition is clear, a mechanism 

^. , of stabilization of the superconducting phase is not apparent. In this paper, we study the dynamical stability of 

the superconducting phase in the gravity dual model of the superconductor and reveal that the superconductor is 

stabilized through a backreaction of the gauge field. 

From the gravity point of view, the normal phase corresponds to charged black holes in anti-de Sitter spacetime, 
i.e., the Reissner-Nortstrom-AdS black holes. The instability below the critical temperature drives the black hole into 
hairy black hole which corresponds to the superconducting phase. Hence, what we would like to prove is the stability of 
this hairy black hole. Therefore, we can use the standard technique in the black hole perturbation theory jU>[l2J . We 
classify general perturbations around the hairy black hole into vector and scalar sectors by means of a 2-dimensional 
rotational symmetry of the black hole. We show the vector sector contains no unstable mode. In order to complete 
the proof of stability, we examine the scalar sector. By looking at the vicinity of the transition point, we demonstrate 
how the stability is realized in the superconducting phase. We also argue that the stability of the system persists 
away from the critical point. 

The organization of this paper is as follows. In section 2, we introduce the model and present background equations 
of motion. Here, we describe the phase transition. In section 3, we explain the Arnowitt-Deser-Misner (ADM) 
formalism which is useful to perform efficient calculations and the stability analysis. In section 4, we prove the 
stability of the vector sector by constructing the positive definite Hamiltonian. In section 5, we reveal a mechanism 
for the stabilization of the superconducting phase through the analysis of the scalar sector. Although the actual 
analysis is performed in the vicinity of the critical point, we argue the stability holds even away from the critical 
point. The final section is devoted to the conclusion. 



II. GRAVITY DUAL OF SUPERCONDUCTORS 



In this section, we review the gravity/superconductor correspondence [3j,[5|. Here, we take into account backreaction 
for completeness. 

The action is given by 
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where 1Z is the 4-dimcnsional Ricci scalar and L denotes the AdS curvature scale. Here, we have incorporated the 
charged scalar field ip and the gauge field A^ from which we can calculate the field strength F^ = d^Ay — d u A^. As 
for the potential V(ip), we take the mass term, V(ip) = m 2 \ip\ 2 . In this paper, we assume that the mass is always 
above the BF bound. Note that the coupling constant q controls the strength of the backreaction. 
Let us consider the static background. Then, the metric is given by: 

ds 2 = -f(r)e-^dt 2 + Tr ^ + - 1 5 ab dy a dy b , 
f(r) L 2 

where a, b denote x, y coordinates. The other fields are expressed by: 
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ip = ip( r ) 
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Note that the scalar field is taken to be real using the U(l) gauge transformation. Then, tt and rr components of the 
background Einstein equations yield 
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Here, a prime denotes derivative with respect to r. while the gauge and the scalar equations become 
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From the last equation, we see the gauge potential (f> induces the effective negative mass squared for the charged scalar 
field. This acts as a trigger of the phase transition. 

By solving the above equations with appropriate boundary conditions, we obtain the asymptotic behavior 
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where /i and p are interpreted as a chemical potential and charge density of the dual theory on the boundary. (Oa) 
represents the expectation value of the operator Oa dual to the charged scalar field ip. The exponent A is determined 
by the mass as A = 3/2 + -\/9 + Am 2 L 2 /2. This is the dictionary of AdS/CFT correspondence. When Hawking 
temperature is above a critical temperature, T > T c , the solution is given by the Reissner-Nortstrom-AdS black holes 
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where r + represents the horizon radius. For T < T c , the Reissner-Nordstrom AdS black hole solutions become unstable 
and new hairy black holes will be firmed. There, the expectation value (Oa) has a non-trivial value proportional to 
\JT C — T . This is the phase transition from a normal phase to a superconducting phase. 

In subsequent sections, we will investigate the stability of the hairy black holes corresponding to the superconducting 
phase. 



III. USES OF ADM FORMALISM 

Although there are some indirect evidences of the stability of holographic superconductors [ljl 13 ; ^ is important to 
give a direct evidence. To this end, we need to perform perturbativc analysis of the system introduced in the previous 
section. In this section, we explain a useful method to obtain a quadratic action for black hole perturbations |15| . 

Since the background spacetime is static, it is useful to work in the ADM formalism. Let us take the parameterization 
for the metric 



ds 2 = -N 2 dt 2 + h l0 (dx l + N i dt)(dx j + N j dt) . 



(10) 



where i = r,x,y denotes the spatial coordinates, N is the lapse function, and N l is the shift vector. Under this metric 
parameterization, the Einstein-Hilbert action can be written in terms of the spatial curvature R and the extrinsic 



curvature Kij as 
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where for convenience we used the symmetric tensor Eij instead of the extrinsic curvature, which is dchncd by 
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Here, a dot and a semicolon are a time derivative and a covariant derivative with respect to hij, respectively. Similarly, 
the action for the gauge field reads 
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The action for the charged scalar field is also written as 
5. = / d^x^g [ - | W - iqA^\ 2 - V(ip) 
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Now, it is easy to see why the ADM formalism is useful for our purpose. In the ADM formalism, the background 
metric reads 



ds 2 b = -N 2 {r)dt 2 + h lJ {r)dx i dx j , 



where 
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Apparently, there is no shift vector N l in the background metric. Hence, we find that E^ vanishes for the background. 
And, the gauge field has only time component, so we have Fy — 0. Thus, some of the terms in Eqs. (|11[) . (fT3")) and 
(1141) are already the second order quantities. In this way, the ADM formalism makes calculations for obtaining 
the quadratic action easier. In order to further simplify the calculations, we take the variation of the total action 
S = Sr + Sa + Sip with respect to the lapse function N to yield the Hamiltonian constraint equation: 
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Substituting this equation into the total action, we can simplify the part proportional to I yhN ) 
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where index (2) means the second order quantity. 



Since there exists a 2-dimensional plane symmetry in the black hole background, we can decompose any tensor into 
vector and scalar sectors by means of the 2-dimensional rotation symmetry. Those sectors are decoupled in the linear 
equations. Hence, we will consider the vector and the scalar sectors separately. 



IV. VECTOR SECTOR 



In this section, we will prove the stability of the vector sector of perturbations. To this aim, we use the formalism 
explained in the previous section. Hereafter, we set L = 1 and k 2 = 1. 
Let us consider the metric perturbations on the background metric (|15l) . 

ds 2 = ds'l + dg^dafdx" . (18) 

The vector sector of metric perturbations is generally expressed by 

/ Sg ta \ 
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V * * Sg ab J 

where Sgta, Sg ra and 5g a b = %C(a\b) satisfy the divergence free condition: <5.g* a ' a = Sg ra } a = Ca'° = 0. Here, | represents 
a partial derivative. Using the gauge transformation x^ — >• x^ — ^ with 



^ = (0,(Ua), 



where £j a = 0, those variables transform into 
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Thus we can make Sg ab vanish by choosing £ a = — £ a and we get the complete gauge fixing. This is called the 
Regge- Wheeler gauge. In the Reggc- Wheeler gauge, the perturbations that belong to the vector sector can be written 
as 

/ v a 

Sg^ v = * w a 
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where vj a = wj a = 0. In the ADM formalism, this corresponds to 
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As to the gauge field, we can take 

5A fi = (0,0, Z a ), (24) 

where Z a satisfies Zj a = 0. 

Now, we can calculate the quadratic action from Eqs. ([TTT) . (fT3")l and (fH)l . The quadratic part of Eq. (jTTJ) is going 
to be 
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The quadratic part of Eq. (fl3[) becomes 
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The total quadratic action for the vector sector becomes 
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where we moved on to the Fourier space with respect to x, y coordinates. Note that there are three unknown variable 
v a ,w a and Z a , which satisfy the transverse conditions. Among them, v a is not a dynamical one which should be 
eliminated. 

The action Eq. (f2"8")l becomes, 
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where we eliminated some terms proportional to w a w a using the background equation @. Now we want to eliminate 
the non-dynamical field v a by completing the square with respect to v a . In doing so, we need the Hamiltonian 
formalism in order to eliminate v' a v' a . Defining the conjugate momentum 
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we obtain the Hamiltonian 
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Let us see the terms containing v a . After integrating by parts for v' a , we have 
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the variable v a can be eliminated because the first squared term vanishes after substituting the equation of motion 
for v a - Then the Hamiltonian becomes 
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We find the Hamiltonian consists of positive terms. This implies the stability of the system. The argument is as 
follows. Suppose an unstable mode exists, then the Hamiltonian must grow. However, since the Lagrangian has a 
time translation invariance, the energy of the system is conserved. This contradicts the growth of the Hamiltonian. 
Therefore, the unstable mode cannot exist in the vector sector. Note that this is true either for the normal phase or 
the superconducting phase. 



V. SCALAR SECTOR 



In the previous section, we have proved the stability of the vector sector for both the normal and superconducting 
phase. Thus, if unstable modes exist, they must belong to the scalar sector. Above the critical temperature, the 
system is the Reissner-Nordstrom-AdS black hole with a trivial scalar field. For this case, we know the system is 
stable [jj| . Below the critical temperature, the effective mass of the charged scalar field violates BF bound, then the 
system ceases to be stable. Eventually, a hairy black hole is formed. It is believed that the new hairy black hole is 
stable. In this section, we reveal the stabilization mechanism of the hairy black hole. 

The scalar sector of metric perturbations is generally expressed by seven degrees of freedom 
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Using the gauge degree of freedom of £t, we can eliminate Sgt a - Using the gauge degree of freedom of £ and £ r , we 
can take Sg a b = 0. Therefore, we take the following gauge. 
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where the components of metric perturbations w,H,Hi and H depend on (t,r,a). In the ADM formalism, we can 
put 
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In addition to the metric perturbations, we need to consider fluctuations of the gauge field SA^ and the charged scalar 
field Stp, Sip*. As to the gauge field, we can take 
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Using the £/(l) gauge transformation, we can take the scalar field perturbations to be real Sip* = Sip. However, for 
the trivial background tp = 0, this gauge fixing is singular because the phase loses its meaning when the amplitude is 
zero. 

Now, we can calculate the quadratic action for perturbed quantities using the ADM formalism. The action for the 
gravity is calculated as 
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In vacuum case, x vanishes and the above action gives the Zerilli equation |ll| . In the present case, we have other 
fields. The action for the gauge field becomes 
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Since we have eliminated the phase of the charged scalar field, the corresponding physical degree is absorbed by the 
gauge field. Then, a and /3 are dynamical degrees describing the massive gauge field. The time component of the 
perturbed vector A is a Lagrange multiplier and hence not dynamical. Finally, the action for the scalar field is given 
by 
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Note that we have eight unknown variables H, Hi, H, w, A, a, f3, Sip. Among these, only four variables are physical. 
Thus, the total action in Fourier space becomes 



s = \ I 



dtdrd 2 kr 2 e 2 



X -r 2 e-^H U^-H' + %H + 2^H + 2^-H - 2^fw' - ^fw - 2^fw 
rf (I f 3 A x , k 2 J ( f 2 X x TT 

k 2 
h^ei (w - H^ 2 + 2re^H 1 H + rfx'e^H 2 

2 U2 



l -U-X' + ^(H-H)] +^'-(/3-A)+^ 7 (/3-A) 2 -||e-M/3'- a ) 



dtdrd 2 k- 2 



re 2 



(Stp - fi''Hi\ - 2q 2 e x (fyil) 2 Hia + q 2 <p 2 5ip 2 + q 2 i\?\ 2 + 4q 2 (PiPSip\ - q 2 fip 2 a 2 



^§V/3 2 - ]H a - %H 2 - W 2 H 2 + 2f4>'6i>'H + 2^-^' w5^ - m 2 5^ 2 



,2rj,2 



+ / dtdrd 2 k-r 2 e-% (H-H)[ 2q 2 ^ 2 \ + 2q 2 cf> 2 ip5i/j + ftp' 2 H - 2 ftp' Sip' - 2m 2 ip5ip 



# 2( 7 j.2.1.2 



dtdrd kr e 2 -^-cp W 



-H 2 - -H 2 - -HH 

8 8 4 



(45) 



where Eq. (j4]) is used. 

Now, it is convenient to use the Hamiltonian formalism for getting an insight. Defining the canonical conjugate 
momentum: 
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we can perform the Legendre transformation and obtain the Hamiltonian density 
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r 2 
-2r 2 feri^'HS^' - 2k 2 Je-iip'w8i/j e"* (H - H) {2q 2 § 2 i\)8i\) - 2fi/)'6ip' - 2m 2 ip6ip) 

+r 2 fe~i5ip /2 + k 2 a e~i8ip 2 + r 2 e^^m 2 8ip 2 

-gVe" V*V> 2 - gVe" V [a 2 + 4^<tyA + <f>(H-&)\\ . (50) 

First, we consider the normal phase ip = 0. In this case, we find the scalar field perturbation decouples from the 
other parts in Eq. ((50)) and reads 



U H = ^e-$i^ + r 2 fe-i54>' 2 + k 2 a e~^H 2 + r 2 e-^m 2 5^ 2 - gVe"*^ 2 ■ (51) 

As the Reissner-Nordstrom-AdS black hole is known to be stable, the instability of the system should be caused by 
the scalar field. Apparently, the last term in Eq. (|5i"j) could destabilize the system. Indeed, as we lower the Hawking 
temperature, the potential <ft oc \/r becomes larger near the horizon. Then, the effective mass m 2 ff = m 2 — q 2 <fi 2 gets 
smaller. Eventually, below the critical temperature, the effective mass m 2 ff = m 2 — q 2 <p 2 violates the BF bound and 
the system becomes unstable. 

Now, we see what happens in the superconducting phase ip ^ 0. In order to make the discussion clear, we focus on 
the stability of the region near the critical point where ip is close to zero and the Reissncr- Nordstrom- AdS black hole 
becomes a good approximation. The scalar perturbation Sip decouples from the metric perturbation in the Rcissncr- 
Nordstrom-AdS limit ip — > 0. Hence, in the limit, we can set H = Hi=H = w = when we look at the scalar 
field. In a sense, this is close to the probe limit. Note that we have four variables a, (3, A, Sip among which only A is 
unphysical. Now, the relevant part of the Hamiltonian is given by 



H = / drdrk, 



i e_iP " + 2*f~* *% + i e_fP ^ + ¥ e ~* (/? ' " ^ + ^/«" Vo» + k 2 a q 2 e-H 2 p 2 
+r 2 fe~i5ijj' 2 + k 2 a e~^Sip 2 + r 2 e~im 2 Sip 2 



-gVe"* W - gVe" V ( A 2 + 4<tyA + f° fi A 

W q 1 r l e ~ip 2 



(52) 



At first sight, the system seems to be unstable because of the term, —q 2 r 2 e~^(p 2 5ip 2 in the last line of Eq. (|52|) in 
the limit of the Rcissncr-Nordstrom AdS ip — > 0. However, ip is not actually zero. We can eliminate the unphysical 
variable A by completing the square for A or alternatively by using the equation of motion for A 

A = -2>- 9 f°r? /2 - (53) 

ip 2q 2 r 2 e ^ip 2 
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Then, the last line of Eq. ([52]) becomes 



+ 3 g Ve- V<V 2 + 2^ {P' a P P ) H + [ P t 2 P ^ 2 ■ (54) 

The effect we have incorporated is merely the backreaction of the gauge field. Notice that the constraint equation, 
which is derived by the variation with respect to A, 

PL - Pp = (55) 

holds in the limit ip — > 0. This constraint equation is nothing but the Gauss law for the gauge field. This appears 
because the limit ip — > recovers the gauge invariance in the gauge field system. Therefore, in the vicinity of the 
critical point, we can ignore last two terms in Eq. (|54[) . Thus, when we take the limit tp — > from the side of the 
superconducting phase, the Hamiltonian reduces to 

H = l^fe^-Le-lp^-^e-fPl + ^e-iP^ + M/e-i^'-a) 2 

+r 2 fe-iSip' 2 + kle'iSiP 2 + r 2 e -*m 2 ^ 2 + SgVe"*^ 2 ^ 2 1 . (56) 

Now, it is clear that the Hamiltonian is positive definite except for the mass term. Moreover, the effective mass 
squared m 2 ff = m 2 + 3q 2 cf> 2 is always above the BF bound. Thus, we find the system is going to be stable after the 
phase transition. It should be stressed that the flip of the sign of the term proportional to <j) 2 5ip 2 is possible only in 
the condensed phase with tjj ^ everywhere. In fact, there is an (infinite) set of solutions for tp which are regular 
on the horizon and satisfy the required boundary condition at infinity. Based on our analysis, we believe that the 
solutions with nodes are likely to be unstable and only the lowest solution with no nodes is stable. 

We can argue the stability of the superconducting phase far from the critical point as follows. From the point of 
view of the gravity dual, the phase transition is the bifurcation of the solutions at the phase transition point from 
where two branches are developed. The normal branch is apparently unstable. What we have shown is the stability of 
the other branch in the vicinity of the bifurcation point. However, from the continuity, this is sufficient to prove the 
stability of the superconducting phase provided that there exists no other bifurcation point on this branch. Indeed, it 
is difficult to imagine a new bifurcation point because the sign flip of the term proportional to cf> 2 Sip 2 persists as long 
as the scalar field has the expectation value. Admittedly, this argument does not give a proof of stability far away 
from the critical point. 

VI. CONCLUSION 

We have studied the dynamical stability of superconductors using the ADM formalism. First, we proved the stability 
of the vector sector by explicitly constructing the positive definite conserved Hamiltonian. We have also obtained 
the Hamiltonian for the scalar sector. In the normal phase, the action decoupled into two parts. One is that of the 
Rcissncr-Nordstrom-AdS black hole and the other is that of the charged scalar field in the Reissner-Nordstrom-AdS 
background. It is well known that the former system is stable. The stability of the latter system depends on the 
effective mass of the scalar field. Since the effective mass becomes more negative as the gauge potential becomes 
large, the phase transition occurs at the critical temperature. And, below the critical temperature, the Reissner- 
Nordstrom-AdS black holes are unstable. Numerically, we know the scalar hair is developed below T c . For this case, 
the structure of the constraint for A changes. After eliminating the unphysical variable A, we found the effective mass 
became above the BF bound. In particular, near the critical temperature ip — > 0, we have decoupled equations for the 
charged scalar field, which is now stable. Thus, we have revealed a mechanism for the stabilization of superconducting 
phase in the scalar sector. The result is remarkable because the stability near the critical point is guaranteed by 
the non-gravitational effect. In fact, the effect of the gauge field changes the signature of the effective additional 
mass when the charged scalar field has condensation. In our argument, we have assumed the solution for the scalar 
field has no node. Therefore, other solutions with nodes are expected to be unstable. Admittedly, we have not been 
able to prove the stability of superconductors far away from the critical point. However, the fact that the effective 
mass squared remains to be above the BF bound in the superconducting phase suggests the stability of holographic 
superconductors. It should be stressed that the stabilization mechanism we found originates from the structure of the 
minimal gauge coupling, hence the result applies to holographic superconductors in any dimensions. 

Our analysis can be extended to various holographic superconductor models. For example, it is possible to extend 
our analysis to the Gauss-Bonnet superconductors J17H21J . In that case, again, we can use the results in black hole 
perturbations [22|, [23| . 
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More importantly, we need to study the dynamical stability of holographic superconductors embedded into string 
theory [24l426l |. However, it requires a more sophisticated method for the analysis. 
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